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Inclusive Quark Production in e + e -Annihilation - 
A Path Integral Approach 

O. Nachtmann, A. Rauscher a * 

a Institut fur Theoretische Physik der Universitat Heidelberg, Philosophenweg 16, D-69120 Heidelberg 

The single-particle inclusive differential cross-section for a reaction a + b — > c + X is written as imaginary part 
of a correlation function in a forward scattering amplitude for a + b ^ a + b in a modified effective theory. In this 
modified theory the interaction Hamiltonian Hi equals Hi of the original theory up to a certain time. Then there 
is a sign change and Hi becomes nonlocal. This is worked out in detail for scalar field models and for QED plus 
the abelian gluon model. A suitable path integral for direct calculations of inclusive cross sections is presented. 



1. INTRODUCTION 

Inclusive cross sections like 

e+ + e" -> n ± + X 

are measured at a centre-of-mass-energy y / i=91.2 
GeV for example by the DELPHI detector at 
LEPQ. So far, only phenomenological models, 
which need to be tuned to the data, are available 
to describe the fragmentation of quarks and glu- 
ons into hadrons. A direct QCD calculation for 
these inclusive cross sections seems difficult be- 
cause of the importance of non-perturbative ef- 
fects. In order to be able to include these soft 
processes, one should start from an expression 
for the inclusive cross section which is suitable 
for non-perturbative calculations. 

The article is organised as follows. In section 
2 the basic relations for single inclusive cross sec- 
tions are recalled. In section 3 we present our 
general formalism for the modified effective the- 
ory in the case of scalar fields. As a specific ex- 
ample the cross section e + + e~ — > q + X is con- 
sidered in quantum electrodynamics (QED) with 
massive photons and the theory of quarks inter- 
acting with abelian gluons in section 4. A path 



* supported by the German Bundesministerium fur Bil- 
dung und Forschung (BMBF), Contract no. 05 7HD 
91 P(0), by Studienstiftung des Deutschen Volkes and 
by Deutsche Forschungsgemeinschaft (DFG), Grant no. 
GKR 216/1-9 



integral representation for inclusive cross sections 
in the abelian gluon model is derived. 

2. SINGLE INCLUSIVE CROSS SEC- 
TIONS 

In this section we recall some basic relations 
for inclusive cross sections. Our notation follows 
HQ. Let us consider a single-particle inclusive 
reaction, i.e 



a(pi) +b( P2 ) ^ c( P3 ) + X( Px ). 



(1) 



To take a simple case, let a, b, c be spinless par- 
ticles with masses m a ,mb,m c . The cm. energy 
squared is s = (pi + P2) 2 - The S- matrix element 
for the reaction (Q) is given as 

s fi = ( c(p 3 ),X(p x ),out I a(px),b(p 2 ),in ) 
= iZ- 1 ! 2 I d A x e txp3 (CL. + m 2 } 



X(p x ),out | <j) c (x) | a(p 1 ),b(p 2 ),in 



(2) 



Here we have applied the reduction formula for 
particle c in the final state. Let <j) c {x) be a suit- 
able interpolating field for c and Z c the corre- 
sponding wave function renormalisation constant. 
The T-matrix element is obtained from the S- 
matrix element via 

Sfi = Sfi + i(27r) V 4) (pi + p 2 - p 3 - Px)Tfi 



Tfi = Z- 1 ' 2 (a y + ml) ( X(p x ),out \ 

Mv) I a(pi),b(p 2 ),in)i . (3) 

The single-particle inclusive cross section f inc (p 3 ) 
is defined by 

d 3 a 
finciPz) := p§-^(a + &^c + X) 

= i y 

A(2TT) 3 w(s,m 2 a ,ml) ^ 

(27r) i 6^(p 1 +p 2 -p 3 -p x )\T fl \ 2 . 

(4) 

In the usual way the sum over all states \X, out ) 
in (0) can be carried out using completeness and 
translational invariancc 



finc(P3) = 



1 



4(2tt) 3 w 



a c 



Z7 1 / d 4 a; e V3a; 



(□« + m=) 



(a(p 1 ),b(p 2 ),in \(f}\{y)(f) c {x + z) 
\a(pi),b{p 2 ),in) ]y ^_^ o _. 

(5) 
Rewrite fi„ c (P3) to get time-ordered operators 



finc(j>3) = 



2(2ir) 3 w 



Im C{p 1 ,p 2 ,p 3 ), 



(6) 



C(px, P 2,p 3 )=i d A xe^ x M(x) 



(7) 



M(s) = ^(□v + m^^ + m?)^-* ) 

(a(p 1 ),b(p 2 ),in \ </>t(y)<f> c (x + z) 
\a{p 1 ),b(p 2 ),in ), 

(8) 

The amplitude C(pi } p2>P3) w ^ be written as a 
field-field correlation function in a forward scat- 
tering amplitude, but in a modified effective the- 
ory. 

3. MODIFIED EFFECTIVE THEORY 
FOR SCALAR FIELDS 

Let us assume that the basic dynamical vari- 
ables of the original theory are the operators for 



unrenormalised scalar fields <j>i{x) and their con- 
jugate canonical momenta IL(x) (i — 1, ...,N). 
We denote <^(x),IL(a:) collectively as $(x). Let 
H be the Hamiltonian of the system which we 
split into a free part Hg and an interaction part 
Hi which may depend explicitly on the time t, 
but should not involve time derivatives of IL (x) 

H{t, $(£ , t)) = H ($(2, t)) + Hi{t, *(£, t)) (9) 

with 

lim H(t,$(x,t))=H ($(x',t)). (10) 

t — >±oo 

Besides the interacting fields and momenta $ free 
fields and momenta <J>(°) are considered with the 
corresponding Hamiltonian Hq- 

We assume now as usual that there exist uni- 
tary operators U(t) that realize the time-depen- 
dent canonical transformations relating $ to $^ - ) 

<S>{x,t) = U- 1 {t)& Q Xx,t)U{t). (11) 

Taking as boundary condition 

lim $($,t) = $ (0) (£,£) (12) 



we get 

U(t) = Texp 



*'#/(£', $ (0) (£,£')) 



We define furthermore 
U(t 2 ,t 1 ) = U(t 2 )U- 1 (t 1 ). 
so that the 5*-matrix is given as 
S = lim U(t,-i) 

t— > — OC 



(13) 



(14) 



= Tcxp 



dt'JTj(t', $ (0) (£,£')) 



(15) 



Now we return to the single inclusive cross sec- 
tion, where we have to calculate the matrix ele- 
ment 



M(x) = Z- c 



□, 



(°« 



9(-x°) 



(a(pi),b(p 2 ),in \ 4>\{y)(j) c (x + z) 
\a(pi),b(p 2 ),in 



la->o-, z^o- 



(16) 



Following the time-dependence in M.(x) from the 
right to the left, we start at time — T — > — oo 
and pass through operators of increasing time ar- 
guments until time 0. Then the time sequence 
changes and we go back in time to time — T" — v 
— oo. In a usual matrix element the time argu- 
ments should increase instead. We will now show 
that we can write the matrix element M(x) in 
the usual form, with time increasing from right 
to left, if we pass to operators U of the form (13) 
but with a modified interaction Hamiltonian Hj . 
As the time-sequence in M. is correct up to t = 0, 
we request 



(17) 



1) Hi(t) = H/(i,$ (0) (f,t)) for t< 

2) 5=1 


with 

S = T exp 

It follows 


r+oo 

-i dt'Hi(t',& 0) (x,t')) 

J — oo 


Hj(t) = 





(18) 



(19) 

For t > our modified interaction Hamiltonian 
Hj (t) depends on the free fields and momenta at 
time (— t). But we know how to express the free 
fields and momenta at time (— t) by their values 
at time t using the free field equations of motion. 
Thus, we can consider Hi(t) as a nonlocal func- 
tional of the dynamical variables <$>(°'(£, t) at the 
same time t also for t > 



H I (t) = H I (t,$W(Z,t)). 



(20) 



Therefore, we can consider our matrix element 
as of the standard type but in the modified theory 
governed by the total Hamiltonian 

H(t, *(g , *)) = tfo($(£ , *)) + #/(*, *(f , *)), (21) 

Here we denote by (( )) matrix elements in the 
modified theory 

M(x) = Z- 1 9{-x°)(a y + m 2 c )(a z +m 2 c ) 
(( a(p 1 ),b(p 2 ),out | </)l(y)(t) c (x + z) 
I a{pi),b(p 2 ),in )), 



The matrix element M. is written in the stan- 
dard form with an in-state to the right and an 
out-state to the left. The prize we have to pay 
is that we have to use the modified theory where 
the Hamiltonian H(t) has a sudden variation at 
t = and is nonlocal for t > 0. 

4. INCLUSIVE PRODUCTION IN QED 
AND AN ABELIAN GLUON MODEL 

4.1. QED and the Abelian Gluon Model 

In this section inclusive reactions are consid- 
ered in QED coupled to an abelian gluon model. 
An example for such a reaction is 

e+( Pl ) + e-(p 2 ) - q( P3 )+X. (23) 

Starting point is the Lagrangian describing the 
interaction of electrons of mass m and charge — e 
with a massive photon of mass A - to avoid any 
infrared divergences - and of two quark flavours of 
equal mass M and electric charge eQ q with the 
photon and a massive abelian gluon of mass \x. 
As Lagrangian we choose 






vlG^ 



KA^ 



V> 



-m + e 4 ip 



+q I - P -M - e Q 9 4- g0 T 3 (?J q, (24) 

where G^ denotes the abelian gluon field and 
G>„ = dfj,G u — d v Gfj, its field strength tensor, 
A^ the photon field and F^ v — d^A v — d v A^ L its 
field strength tensor, ip the electron field and q 
the quark field. We have a quark field with two 
flavours and T3 is the usual Pauli matrix. 

The Hamiltonian corresponding to ( g4|) can be 
obtained as usual by a Legendre transformation. 
Only the interaction term has to be altered in 
order to construct the Hamilton operator for the 
modified theory 

Hi(y°) = fd 3 y [Oi-y^Hjiy) -0(|/°)Wj(»)] 



lj/->0-, z^O- 



with y = (-y°,y). 



(22) 



(25) 



In Hi(y) the fields at time — y° have to be substi- 
tuted by the fields at time y°. Two points should 
be stressed: First, due to the jump in the interac- 
tion, the effective theory has no time translation 
invariance, so that no longer energy conservation 
holds at every vertex. Second, the effective inter- 
action contains derivatives. 

4.2. Inclusive Quark Production in e + e - 
Annihilation 

For unpolarised e~ and e + in the initial state 
and summation over spin and flavour of the final 
state quark the matrix element reads 

M(x) = -Z-^uipa^-M) 

s q 
]r'(( e + ( Pl ),e-(p 2 ) \Tq A (x + z) 

S 4. ,5 _ 

q A (y)\e + (pi),e-(p 2 )))9(-x°) 
Hft-M)u(p3)| 

(26) 

In order to calculate M. we expand it in pow- 
ers of the electromagnetic coupling constant e. 
The electromagnetic interaction of the incoming 
fermions can be separated 

2 

C a (p u p 2 ,p 3 ) = Z^e^^F" 



1 



A 2 + ie 

^u(p 3 )(iK-M)((0\Tq A (x + z) 

!d 4 x 2 [6(-x 2 )e-^+P^^q(x 2 ) 1 ,q(x 2 ) 
-e{xl)e-^ + ^-q{x 2 ) lv q{x 2 )] 

Id^iei-xly^+p^qix^qix!) 

-Oixly^+^qix^qix!)} 
q A (y)\0))(-i% - M)u{p 3 \ 

I y— »0~ ,z^0~ .e=0 

(27) 

As a simple check of our theoretical manipula- 
tions let us finally calculate the lowest order in the 
quark-gluon-coupling g, i.e. for g = 0. Decom- 
posing the quark 6-point-functions with Wick's 



theorem into quark 2-point-functions and substi- 
tuting for these the perturbative propagators we 
get for C a two contributions according to the two 
different possible contractions shown in figure 1. 




Figure 1. Lowest order contributions. 



The distinction between both contributions is 
given by the fact, that the upper diagram has 
got an imaginary part, whereas the lower one is 
completely real so that we get only a contribution 
from the upper diagram 



fincipz) 



1 



2(2tt) 3 u; 



we 4 Q 2 q F v 



1 



A 2 



le 



2j2^(P3hA^-Mh„u{p 3 ) 

s q 

S (1 Hp 2 - M 2 )9(p°) 



Pi +P2 ~P3 

(28) 

which is, of course, the standard result, which one 
obtains in considering to lowest order the reaction 



e7 + e — > q + q 

for two quark flavours of charge Q q . 



(29) 



4.3. Path Integral Representation 

In the last subsection M. and therefore the one- 
quark-inclusive cross-section could be expressed 
in terms of quark 6-point-functions after separat- 
ing the electromagnetic interaction by a perturba- 
tive calculation. As the coupling g is not assumed 
to be small we will now derive a representation 
of these quark 6-point-functions suitable for non- 
pcrturbative calculations. For this we consider 
the Hamiltonian path integral in our effective the- 
ory obtained from the abelian gluon model. 

Any Green's function of the theory can be writ- 
ten as 



((0| Tq(x 1 )...q(x 2 ).. 



(30) 



V(G,Il G ,q,q) q(x 1 )...q(x 2 )... 

xexp{iJd\(Il(y)<i>(y)~H(y))} 
with 
Z = f-D(G,U G ,q,q) 

x exp{ ? JdSj (n(y)^(y) - H(y))}. (31) 

Here the classical fields G and H G and the Grass- 
mann fields q and q have to be inserted into the 
part of 7Y(y) d2q) which describes the abelian 
gluon model. Furthermore for y° > the classical 
and Grassmann fields have to be time-shifted. 

As the Hamiltonian path integral is quadratic 
in q, q, the fermionic fields can be integrated out. 
With the Green's function Sf of the quark in the 
modified gluon background 

(i^M Q -g&(z 1 ))S F (z 1 ,z 2 ;G) 

= -S^( Zl -z 2 ), (32) 

the quark propagator is given as 

<( | Tq( Zl )q(z 2 ) | 0)) = « -S F ( Zl , z 2l G))>, (33) 

i 

where the brackets (( )) on the right hand side 
denote the average over all gluon fields with the 
measure dictated by the path-integral 

((F(G, n G )» = i- 1 J v(G, u g )f(g, n G ) 

x cxp{i (dSj (n(i/)#(i/) - H(y))} (34) 



Inserting the Green's function into the matrix 
element we get the inclusive cross section ex- 
pressed explicitely in terms of a Hamiltonian path 
integral. This expression should be a conve- 
nient starting point for applying non-perturbative 
methods to an evaluation of inclusive cross sec- 
tions. It should be possible, for instance, to gen- 
eralise the methods of || to the case of this effec- 
tive theory here. 

5. CONCLUSIONS 

We have written the inclusive cross section 
as imaginary part of an amplitude C for which 
we have given a path integral representation. 
We have applied our formalism to the reaction 
e + + e~ — > q + X and checked that in lowest non- 
trivial order we get the correct result. Our hope 
is that our path integral representation will lead 
to a genuinely non-perturbative evaluation of in- 
clusive cross sections at high energies along the 
lines of Q. The generalisation to QCD should 
be straightforward. Our methods should allow a 
description of inclusive production of hadrons h 
in e + e _ annihilation at high energies 

e + + e - -»• h + X (35) 

in terms of a genuine OPE. A similar approach 
should be possible for fracture functions M for 
hadron-hadron scattering 



hi + h 2 



h 3 +X. 



(36) 
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